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Abstract - A tilting phase transition is predicted for systems comprising 
rod hke molecules which are irreversibly grafted to a flat surface, so that the 
non interacting rods are perpendicularly oriented. The transition is controlled 
by the grafting density p. It occurs as p increases as a result of the interplay 
between two energies. Tilt is favoured by the van-der-Waals attraction between 
the rods. It is opposed by the bending elasticity of the grafting functionality. 
The role of temperature is discussed, and the tilting mechanism is compared 
to other tilting transitions reported in the literature. 
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1 Introduction 

Anisotropic molecules forming thermotropic and lyotropic liquid crystals can 
exhibit tilt. The onset of tilt is the characteristic of the smectic A-smectic 
C phase transition []T|. Tilted phases are observed in monolayers formed by 
amphiphilic molecules^]. A variety of tilting phase transitions are possible at 
solid-fluid interfaces incorporating anisotropic molecules Diverse molec- 
ular mechanisms can give rise to such phase transitions. In the following, we 
propose a novel mechanism for the development of tilt at surfaces carrying 
grafted rods. In particular, we consider a fiat solid surface supporting cova- 
lently bound rodhke molecules at a surface density p < 1. The rods are grafted 
to the surface so that isolated rods are oriented along the surface normal. The 
grafting functionality, the "joint", is assumed to allow bending at the price of 
an elastic energy penalty. A tilting phase transition is predicted to occur as the 
grafting density p increases. The tilt is favoured by the van-der-Waals attrac- 
tion between the rods, but is opposed by the bending penalty. The interplay 
between these two energies gives rise to a tilting phase transition. The transi- 
tion is predicted at zero temperature. However, the mechanism still works at 
sufficiently low finite temperatures, as we shall argue in the discussion. 

From a fundamental perspective, this mechanism is of interest as a com- 
plement to the variety of tilting mechanisms studied in the literature. Tilting 
transitions are also of importance vor the design of alignment layers for liq- 
uid crystal displays[§]. The orientation of the nematic director inside a liquid 
crystal cell with no electric or magnetic field applied is usually determined by 
the interface. Properly designed interfaces supply the desired in-plane asym- 
metry as well as a certain "pretilt" angle. Grafted liquid crystal (LC) layers 
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may form good alignment layers, when the LC-director adopts the tilt of the 
grafted layer. In practice, the tilt of the nematic director may not be equal 
to the tilt of the grafted layer due to the complicated interaction between the 
nematic medium and the layer. This does not however detract from the prac- 
tical potential of this system. A tilting phase transition driven by the grafting 
density may provide a scheme to fabricate alignment layers with an arbitrary 
tilt angle. 

Another aspect of grafted layers of interest to liquid crystal alignment is 
their behavior under nematic orientational stress. Such stresses occur during 
the Frederiks transition and are of high practical importance as well. Coupling 
between the tilting phase transition and the Frederiks transition is expected 
to reduce the switching voltage of such devices^]. 

There are a number of reports, where rods have been covalently attached 
to a surface |T^, 0]. In most cases, however, the joint is assumed to be 
so flexible that it will not show any resistance to bending. Recently, attach- 
ment with rigid joints has been claimed as well[ll2|, |13|. Note however, that 
a rigid joint does not necessarily imply perpendicular orientation. The angle 
of the joint will depend on details of the chemistry. Right-angled joints are 
the exception rather than the rule. If, however, the tilt angle is small and the 
in-plane direction of tilt is random, it is reasonable to describe the rods as 
perpendicularly oriented. 

The essential physics involved is discussed in section II for a highly simplified 
case: A linear array of grafted rods. An analysis of the two dimensional case 
is presented in section III. A comparison to previously discussed tilting phase 
transitions is given in the discussion. 
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2 The Linear Array 



The underlying physics of the tilting transition are most easily understood in 
the case of a linear array. The rods, of length L and diameter a, are grafted 
onto a straight line at regular intervals of width D. In the absence of rod-rod 
interactions, the grafting functionality imposes perpendicular orientation of 
the rods, i.e. the angle between the rods and the surface normal is 7 = 0. 

However, the joint can be deformed with a bending elastic penalty of 

00 

f/feend = 51 ^isin2^7, (1) 

where Ai are non negative phenomenological constants specified in units of 
energy. The rods also experience mutual van-der-Waals attraction. We picture 
the rods as comprising monomers of size a|T^ . The non retarded van-der-Waals 



attraction between two monomers separated by a distance r is 

w = -c/r^, (2) 

where c is a material constant with units of [energy] [length]^. The van-der- 
Waals attraction between two parallel rods, E{D,'y), may be obtained by 
integration of the pairwise interactions between monomers belonging to the 
different rods (see Appendix A). At D/L < 1, the leading terms in E{D,'~f) 
are 

E{D,j) = - ^ + t/cos7(cos2 7 - 3 - 37tan7) + 0{d')} (3) 

(c(cos7j'' 2 

with d = D/L and K = c/{Aa'^L^). Note that E{D,'~f) favors large tilts since 
in this limit the reduced distance between the rods is Dcos'j < D. Notice 
further that E{D,'~f) diverges at 7 = 7r/2, when the rods are fully tilted. This 
last feature is however an artefact, since 7 = 7r/2 is actually unattainable 



because of excluded volume interactions between the rods. An upper bound 
for the maximal possible tilt in the D/L < 1 case is 7„p = arccos(a/D). 

To analyze the phase behaviour, it is convenient to consider E{D,'y) in 
the limit of ^ 1 and 7 -C 1. In this limit, where corrections due to the 
finite length of the rods are negligible, the leading terms in the expansion of 
of E{D,'j) in powers of sin 7^ are 
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E{D,'j) = -{37TK/2d^)[l + (5/2)sin27 + (35/8)sinS + ■ ■ ■] = - ^ S^sin^S- 

i=l 

(4) 

We now consider the total van-der-Waals energy, U^dw, of a rod in a uniformly 
tilted array. To obtain U^dw we sum the pairwise interactions of the rod i.e. 
Uvdw = 2 E^=i E{nD, 7). In the limit of D/L < 1, when all terms in E{D, 7) 
exhibit a dependence, the result is 

U,dw = 2C{5)E{D,^), (5) 

where ({5) = ~ 1.037 is the appropriate Riemann zeta function. In 

doing so, we assume that the van der Waals interactions are additive. That 
is the case when the dielectric susceptibility of the materials is an additive 
function of the constituting components i.e. when the susceptibility is pro- 



portional to the concentration of the rods[15]. One should however note that 



this requirement in not always fulfilled [|16[. If one considers as an approxi- 
mation only nearest neighbour interactions, the calculated interaction energy 
Uvdw = '2E{D, 7) is very close to (j^). Hence the total van der Waals energy is 
dominated by the nearest neighbour contributions. 

The behaviour of the linear array of rods is determined by the total energy 
per rod U = Ubend + Uydw- It is instructive to consider the behaviour of U in 
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a number of scenarios. If one overlooks the excluded volume effect between 
the rods, U^ciw diverges at 7 = 7r/2. When Ubend is finite at 7 = 7i/2, this 
divergence dominates the behaviour of U, f/(7 = 7r/2) = —00, indicating that 
all the rods are tilted. While U{'~f) may exhibit two minima, U{tx/2) = —00 
always corresponds to the equilibrium state in this scenario. A richer repertoire 
is possible when the bending penalty at 7 = 7r/2 is infinite as well. In that 
case, all terms in the expansion of Ubend have to be retained. Consider, for 
example, the situation when the coefficients of sin^*7 of Ubend and U^dw are 
such that Ai — Bi > for alH > 2. A second order tilting phase transition 
is expected when Ai — Bi = i.e. Ai/K = 157c/Ad^ with d = D/L or 
d = dc = (157rfr/4Ai)^/^. The rods are perpendicular to the surface, 7 = 0, 
when d > de, while for d < de they are tilted at angle of 7 oc y/de — d. A first 
order phase transition is expected at Ai — Bi > 0, if some of the higher order 
terms are negative and large enough such that J2i>2{^i ~ -Bj)sin^*7 < for 
some value of 7 < jup- 

The role of higher order terms of Ubend is certainly important. A discussion 
based on this effect suffers however from two disadvantages. Typically, the 
higher order Aj's are unknown. Also, a chemical design allowing for control 
of these parameters is impractical. To simplify the analysis, we assume that 
the bending energy is Hookean up to a cutoff angle 7„p. For specifity, we 
assume that it is the result of excluded volume interactions and hence given 
by 7up = arccos(a/i5) as noted earlier. Note however 7„p may also be imposed 
by the bending potential. 

It is helpful to consider first U in the limit of d <^ 1 and 7 <^ 1, retaining 
only the first order terms of Uydw and Ubend as specified by (|lD and (^. U is 
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then given by 

U =[Ai- 157vKC{5)/2d^]sm^'y. (6) 

When e = Ai- 157rKC{5)/2d^ > 0, the system may exhibit a first order tilting 
phase transition. We shall discuss this scenario in greater detail shortly. A 
second order transition is expected when e = 0, while if e < the system will 
be tilted. To analyze the first order transition, it is helpful to express U in 
terms of (|l|) and (^. The energy per rod in the limit of (i -C 1 is then 

[/ = Aisin% - a/((i^cos^7), (7) 

where a = 3'n'({5)K and 7 is allowed to vary between and7„p. U is minimal 
at 7 = 0, where U{0) = —a/d^, and at 7„p, where U{jup) = Aism'^'jup — 
a/i^d^cos^'jup)- It is maximal at < ■jmax < Jup, defined by cos'^7max = 
5a/2Aid^. A first order phase transition between an untilted phase and 
a phase tilted at an angle 7„p occurs when f/(0) = U{'^up) or a/d^Ai = 
sin^7„p/(cos-^7„p - 1). 

3 The Planar Case 

To extend our considerations to two dimensional, planar arrays, it is neces- 
sary to allow for two geometrical effects. First, for a given rod the number 
of interacting neighbors increases with the distance. Thus, a shell of radius 
r and thickness dr contains prdr interacting rods, while in the linear array 
the number is always pdr. Second, in the two dimensional case the shortest 
distance between parallel rods, 6, varies with azimuthal angle, 0, between the 
vector joining the grafting sites and the direction of the tilt. When = the 
shortest distance is 5 = rcos7 while for = 7r/2 it is 5 = r. In general, 6 is 
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given by S — r(l — cos^0sin^7)^/^ (Figure 1). In realistic situations, it is also 
necessary to allow for the randomness of the grafting sites. These are unHkely 
to form a regular lattice. It is thus useful to characterize the array of grafting 
sites by a pair distribution function 

g{r) = ph{r). (8) 

The introduction of g{r) also simplifies the mathematical treatment of the 
problem. The van-der-Waals interaction energy per rod in a uniformly tilted 
array is 

Uvdw^P rdr d(f)h{r)V{r,(l),-f), (9) 
Jo Jo 

where V^(r, ^,7) E{r,9) and 9 is defined by sin^ = sin 7 cos 0. It is then 
possible to specify U^dw in terms of the inverse moments of h{r) 

fJ'3 — h{r)dr/r^ l^i—j h{r)dr/r'^, (10) 

Jo Jo 

leading to 

U.MPn) = -^^^s[2FAlA,sin'^)^l - ^]+n{^) (11) 

2 2 2 cos^7 

with the hypergeometric function 2-Fi(l/2, 5/2, sin^ 7) = Jq^"^ d0[l— sin^ 7C0S^ (f>\~^^'^ 

and I = L(3/Li4)/(2/i3). The remainder TZ is of order 1/L^ (sec Appendix B). 

For simplicity, we take h{r) to be a simple step function 

h(r) = e(r - D) (12) 

with the minimal distance between rods D. At D <^ L, the inverse moments Hi 
arc then given by /ij = l/((i — 1)D^~^) and I = L/D. We further assume that 
the minimal distance between rods depends on the grafting density according 
to a proportionality law p — v/D'^. This implies that the rod grafting sites are 
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distributed evenly on the plane, like in the case of the linear array discussed in 
section 2, where the distance between rods is kept fixed. The van-der-Waals 
energy can then again be expanded in powers of sin 7. 

U,dw{pn) = -i^I^Vp'{l+(5/4) sin2 7+(105/64) sin^7+---} = -Y,B,sin^'^ 

(13) 

with K = CT^ / [Aa? . After combining this expression with the bending 
energy per rod (|l]), one can follow the discussion from the previous section: A 
second order transition is expected at Ai = Bi, i.e. at the critical density 

P = Pc = {^?" ■ L-^l'. (14) 
oil 

It may be preempted at some lower density by a first order tilt transition, if the 
cutoff angle 7„p is large enough such that V^dw{Pc, lup) - U^dwiPc 0) + Ubendilup) 
The van der Waals energy in the limit D/L <^ 1 scales like L^^^, whereas 
the bending energy per rod is independent of rod length and grafting density. 
Hence the densities, at which phase transitions occur, scale with the rod length 
L like p, oc 

We note that these considerations only hold under the supposition that the 
rods are distributed evenly on the plane, D oc 1/^/p, as mentioned above. 
If the rods are grafted independently of each other, a Poisson distribution is 
expected. In this case, the assumption D = a with the monomer size a may 



be more appropriate. Eqn (|T3]) then has to be replaced by 



U,dw{p, 7) = + (5/4) sin2 7 + (105/64) sin^7 + ■ ■ ■} (15) 

and densities of phase transitions scale with the rod length as pc 1/-^- A 
second order phase transition is expected at pc = {I6a^ Ai) / (bncL) . However, 
in this case, the validity of a mean field averaging as performed in eqn (0) 



becomes questionable. This is because domains with different tilt order and 
local defects of the tilt order parameter may become important even at zero 
temperature. 

4 Discussion 

The tilting behaviour considered in this paper results from the interplay be- 
tween three factors: The van-der-Waals attraction between the rods, the bend- 
ing elasticity of the joint and the grafting constraints. Clearly no joint is per- 
fectly rigid. The necessary grafting constraints involve two ingredients: (i) The 
bonding of the rods' ends to a fiat surface in a manner imposing perpendicular 
orientation of the non interacting rods, (ii) The absence of lateral mobility. 
The rods start to tilt at the minimal grafting density Pc, where the effect of 
the attractive interaction starts to overcome the bending penalty. Upon fur- 
ther increasing the grafting density, the excluded volume interactions between 
rods gain importance. As the rods are pushed together, they gradually stand 
up and finally undergo another continuous transition to an untilted state at 
a — D, where jup = arccos(a/-D) = 0. This second transition is driven by the 
interplay between attractive and repulsive interaction. Thus one expects tilt 
order in a well-defined density regime (Figure 2). 

We should note that, at zero temperature, there is no entropic opposition 
to the in plane symmetry breaking associated with collective tilt. At finite 
temperatures, however, the entropy comes into play, giving rise to a variety of 
new phenomena: 

First, it is well known that Goldstone excitations will destroy any true long 
range order in two dimensional systems with continuous symmetry, where the 
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ground state is continuously degenerated [0]. In our case, the Goldstone modes 
are long wavelength fluctuations of the in plane tilt direction. At sufliciently 
low temperatures, however, one still has quasi long range order, i.e. the cor- 
relation functions decay algebraically. Unless preempted by a flrst order tran- 
sition, Kosterlitz-Thouless-type unbinding of vortex excitations is expected to 
destroy this order at high temperatures. Note that the energy contained in the 
vortices will approach zero as the tilt goes to zero. Therefore systems with a 
small tilt will undergo a Kosterlitz-Thouless transition at lower temperatures 
than systems with high tilt. This implies that the tilting phase behaviour at 
finite, sufficiently low, temperatures is very similar to the zero temperature 
behaviour - but with two differences: (i) the tilted state now has quasi long 
range order rather than true long range order, and (ii) the second order tilting 
transitions are replaced by Kosterlitz Thouless transitions. 

Second, the entropy being another factor which opposes collective tilt, it 
may assume the role of the bending elasticity at finite temperatures. Indeed, 
the competition of attractive interactions between rods and rod orientational 
entropy gives rise to a scenario of tilting transitions which is very similar to 
the one discussed heref^, |^. The details may be different because the entropy 
will depend on the lateral density due to the packing constraints, whereas the 
bending elasticity is independent of lateral density. 

When lateral mobility is allowed, phase separation will replace the tilting 
phase transition (cf. |0, |1^), because the surface free energy of untilted rods 
is higher than that of tilted rods. Note however that this last scenario can be 
modified by several factors. First, if the rods are given some flexibility, the loss 
of surface free energy per chain upon lowering the chain density is compensated 
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in part by the gain of conformational entropy. As a consequence, phase sepa- 
ration occurs at lower surface densities. Self consistent field calculations show 
that, for sufficiently flexible chains, the areal density of the condensed phase 



may be low enough to favour tilt|]T9[. Second, the phase separated condensed 
phase may support tilt even in systems of perfectly stiff rods if this is ener- 
getically rewarded, e.g. by attractive chain-surface interactions!^, |2y, |2l|, |22 



Third, the distance of closest approach between grafting sites may be larger 
than a, the rod diameter. This last situation is reminiscent of the onset of 
tilt in amphiphilic monolayers, when the distance between "grafting sites" is 



determined by the effective radius of the ionic head group [Q. A similar mech- 
anism [^] for the development of tilt has been proposed for lamellae formed by 
rod-coil diblock copolymers i.e. polymers consisting of a rod like block joined 
to a flexible, coil like block. These can form lamellae in a selective solvent 
which is a precipitant for the rods. In this case the surface free energy of the 
close packed rods opposes the tilt which is favoured by the repulsive interac- 
tions between the overlapping coils. These last tilt mechanisms are however 
less closely related to our system since they do not involve rods grafted to a 
surface with p < 1. The rods are always closely packed. Furthermore, the driv- 
ing force for the tilt comes from units outside the rods (substrate, head groups 
and flexible blocks respectively) and will not operate in a system consisting 
solely of rods. 

As discussed in the introduction, tilting transitions are particularly impor- 
tant, from a technical point of view, in relation to liquid crystal alignment. If 
the surface is in contact with a nematic liquid crystal, the tilt of the rods at 
the surface will influence the director in the bulk of the liquid crystal, and vice 
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versa. Tilt transitions such as the one discussed here will then translate into 
"anchoring transitions" , which involve changes of the direction of alignment in 
the bulk of the liquid crystal. 

In this context, it is instructive to briefly discuss the tilting phase transi- 
tion predicted for liquid crystalline polymers (LCPs) grafted onto a flat surface 
immersed in a good nematic solvent 0. In this case the tilt results from the 
interplay between the repulsive monomer- monomer interactions and the elas- 
ticity of the solvent. This situation involves a number of distinct features. 
First, the effect is expected in the case of long, main chain LCPs supporting a 
number of hairpin defects. Such LCPs adopt a highly anisotropic, ellipsoidal 
form such that their major axis is aligned with the nematic director. Second, 
the effect depends on the anchoring conditions imposed by the grafting surface. 
It is predicted when the anchoring is homogeneous i.e. the nematic director 
at the interface is parallel to the surface. As a result of the coupling between 
the LCPs and the nematic solvent, the orientation of the isolated LCPs is with 
their major axis parallel to the surface. As the grafting density increases, the 
overlap between the LCPs grows. This, in turn, results in an increase of the 
repulsive interactions between the LCPs. These can be weakened by tilting 
the LCPs thus decreasing the monomer volume fraction within the layer. The 
onset of tilt is however opposed by the elasticity of the nematic solvent since 
it results in a distortion of the nematic medium. The interplay of these two 
contributions gives rise to a second order tilting phase transition. Tilt is now 
favoured by repulsive interactions rather than by attraction. The role of the 
joint is assumed by the nematic solvent. The coupling between the LCP and 
the oriented nematic determines the orientation of the isolated LCPs and gives 
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rise to the elastic penalty that opposes tilt. 

As compared to the case studied in the paper and the tilting transitions 
discussed so far, the last scenario is different in that the tilting transition does 
not involve in-plane symmetry breaking. The azimuthal symmetry is already 
broken in the "untilted" state, in which the LCPs lie parallel to the surface. 
This is an effect of the nematic medium, which insures collective alignment 
even at the absence of direct interaction between the molecules. Such tilting 
transitions without symmetry breaking are frequently found in liquid crystal 
systems. 

An experimental example are anchoring transitions induced by trans-cis 
isomerizations of azobenzene molecules, which are covalently attached to the 
surface]^, |10|, |Tl| and immersed in a nematic LC solvent. In these cases, the 



rods are attached to the surface via flexible spacers comprising 4-10 CH2- 
units. UV-irradiation induces a transition from a rod like trans conformation 
to a highly bent cis conformation in the azobenzene molecules. The rods in 
the trans state favour perpendicular (homeotropic) alignment of the liquid 
crystal molecules, whereas the substrate and the molecules in the cis state 
favour parallel (homogeneous) alignment. At a high density of trans units, 
the competition of these two factors leads to tilted anchoring. UV-irradiation 
causes the density of trans units to decrease, thereby inducing a continuous 
transition to an untilted state with parallel anchoring^. 

5 Conclusions 

We have described a tilting transition in grafted layers of rigid rods attached 
to the surface via a chemical bond, which resists bending. The competition 
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between bending elasticity on the one hand and van der Waals attraction on 
the other hand drives the transition. The transition can be of first or second 
order. We find perpendicular alignment for low coverage, and tilted alignment 
for higher densities. Second order transitions will be preempted by Kosterlitz- 
Thouless transitions at nonzero temperatures. We compare the transition with 
a number of related tilting transitions, which have been discussed in the context 
of monolayers and of liquid crystal alignment. In some cases similar scenarios 
are found, where the entropy takes the role of the stiff chemical bond. 
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Appendix A: The one dimensional case 

We consider two parallel rods of length L grafted at distance D from each other, 
and uniformly tilted by an angle 7. The sum of van der Waals interactions 
between monomers from different rods is then given by 

E{D,-i)^-(cla^) dx C dx'- -— . (16) 

^ ^ ' .Jo Jo [(Dcos7)2 + (a:-a:' + Dsin7)2]3 ^ ' 

We can decompose the integral and rewrite E{D, j) as E — — c/ (4a^L) • E{D/L, 7) 
with 

1 — X 



%,7) = 8[^'dx- 



/ dx(e — x)^ ^ ''"^^ "'"^^ 

^0 
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where 5 = dcos'j and e = dsiwj. The integrals can be solved exactly, and one 
obtains 

(ct cos 7)* (1 + ct^)"' — 4a^ sm 7 
, 1 + (i sin 7 ,1 + d sin 7 1 — sin 7 1 — d sin 7 , 



+5 



i -t- a sin 7 i — u sin 7 l — u sin 7,11 

, arctanf — ; ) + ( — ; ) arctanf — ; ) }. 

d cos 7 d cos 7 d cos 7 ct cos 7 -i J 

This function can be expanded in different ways for the cases d < 1 and d > 1: 

(a) d < 1: 

E{d, 7) = 7;^ vA^TJ^ + cos' 7 - 3 - 37 tan 7} + 7^l(d, 7) (18) 

[d cos 7)^ 2d cos 7 

with 7^l(rf, 7) = 2/5 + ^2(4 _ 32/7 ■ cos^ 7) + 0{d^) 

(b) d>l: 

m 7) = ^ + ^(14 - 16 cos^ 7) + Oi^) ^ n2id, 7) (19) 
Appendix B: The planar case 

More generally, let the two rods be grafted on the xy-plane and uniformly 
tilted in the y-direction. The vector connecting their grafting points can be 
parametrized as r = r{cos(f), sin0, 0) and the tilt vector is (sin 7, 0, cos 7). The 
interaction between the two rods is then V{r, 0; 7) = E{r, 9) with an effective 
tilt angle 9 determined by rsin^ = rt (see Figure 2). We now consider an 
assembly of rods, which are randomly distributed on the plane according to 
an isotropic pair distribution function g{r) = ph{r). The total van der Waals 
energy per rod, given by eqn (^^, can then be calculated explicitly for the 
leading terms in the expansion ([18|) . Using the identities 

/ — r = — -^-Ti -. -. 1. sm 7), 
io ^/T^~^^ 2 ^^2' 2' ' 
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with z = sin 7 cos (f), and 



V2 ^ r32;arcsin2; ^ 2:^ + 2 



TT 



cos^^' 



one obtains eqn (|lTD. The order of magnitude of the neglected rest 71 can be 
estimated by 



pc 



r27r rl roc 

71= / d(p { dx X h{xL)TZi{Xj ^) + dx x h(xL)TZ2{x, 9)} 
Jo Jo Jl 

rl roo 

< 2TThmax{J dx X {2 / b + Ax^ + ■■■) + j dx X {A/ x^ + U / x^ + ■ ■ ■)} 



= 0{1), 
hence TZ is of order 
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Figures 



Figure 1: Two parallel rods grafted on the xy-plane 
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Figure 2: A schematic representation of the tilting transitions upon increasing grafting 
density p at low temperatures 
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